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ABSTRACT 

Context. The space mission Kepler provides us with long and uninterrupted photometric time series of red giants. This 
allows us to examine their seismic global properties and to compare these with theoretical predictions. 
Aims. We aim to describe the oscillation power excess observed in red giant oscillation spectra with global seismic 
parameters, and to investigate empirical scaling relations governing these parameters. From these scalings relations, we 
derive new physical properties of red giant oscillations. 

Methods. Various different methods were compared in order to validate the processes and to derive reliable output 
values. For consistency, a single method was then used to determine scaling relations for the relevant global asteroseismic 
parameters: mean mode height, mean height of the background signal superimposed on the oscillation power excess, 
width of the power excess, bolometric amplitude of the radial modes and visibility of non-radial modes. A method 
for deriving oscillation amplitudes is proposed, which relies on the complete identification of the red giant oscillation 
spectrum. 

Results. The comparison of the different methods has shown the important role of the way the background is modelled. 
The convergence reached by the collaborative work enables us to derive significant results concerning the oscillation 
power excess. We obtain several scaling relations, and identify the infiuence of the stellar mass and the evolutionary 
status. The effect of helium burning on the red giant interior structure is confirmed: it yields a strong mass-radius 
relation for clump stars. We find that none of the amplitude scaling relations motivated by physical considerations 
predict the observed mode amplitudes of red giant stars. In parallel, the degree- dependent mode visibility exhibits 
important variations. Both effects seem related to the significant influence of the high mode mass of non-radial mixed 
modes. A family of red giants with very weak dipole modes is identifled, and its properties are analyzed. 
Conclusions. The clear correlation between the power densities of the background signal and of the stellar oscillation 
induces important consequences to be considered for deriving a reliable theoretical relation of the mode amplitude. As 
a by-product of this work, we have verified that red giant asteroseismology delivers new insights for stellar and Galactic 
physics, given the evidence for mass loss at the tip of the red giant branch. 

Key words. Stars: oscillations - Stars: interiors - Stars: evolution - Stars: mass-loss - Stars: late-type - Methods: data 
analysis 



1. Introduction 

The CNES CoRoT m ission (iMichel et al l |2006[ ) and the 
NASA Kepler mission (jBorucki et al.ll20 10') provide us with 
thousands of high-precision photometric light curves of red 
giants a nd have opened a new era in red giant asteroseis- 
mology ()De Ridder et al.|[2009| ). The asteroseismic analysis 
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of these data gives a precise view of pressure modes (p 
modes), corresponding to oscillations propagating essen- 
tially in the large stellar convective envelopes, as well as 
of mixed modes, corresponding to pressure waves coupled 
to gravity waves propagating in the core radiative regions. 
Previous work has mainly dealt with the frequency proper- 
ties of oscillations. In this paper, we explore and quantify 
the oscillation energy. This task is of prime importance for 
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examining the link between convection and oscillations, in- 
cluding aspects such as mode excitation and damping. In 
parallel, examining the power density of the background 
signal will provide useful information on granulation. 

The present work builds on previous studies based 
on CoRoT and Kepler observations, which derived first 
insights on the red giant oscillation spectrum and pro- 
vided empirical scaling relations between the global seis- ^ 
mic para meters ( Dc Ridde r et al. 2009; Hckker ct al. 2009t 
Carrier e t al. 2010; Beddi ng et al.il201oi:lM'osser et al.ir20ia 
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Huber et ah 2010: Mosser et ani2011bD . Very relevant to 



this work, iMathur et al. (201^ have determined the back 
ground signal due to the different time scales of activity and 
granulation, for the sa me set of Kepler red giants consid- 
ered here. Additionallv. iBedding et all (|201lf ) have recently 
shown with Kepler observations that the evolutionary sta- 
tus of red giant can be derived from the measurements of 
the gravity-mode spacing observable in £ = 1 mixed modes 
(|Beck et al.ll2011[ ). With CoRoT data. lMosser et all (|2011aD 
have shown that the evolutionary status has also a clear ef- 
fect on the oscillation power excess. 

In this work, we provide a first complete descrip- 
tion of the oscillation power excess and derive the 
bolometric amplitudes of radial modes. As in previous 
work, we establish scaling relations for global parame- 
ters as a function of frequency. Then, in order to pro- 
vide a direct observational compa rison to theoretical pre- 
dictions (e.g. iDupret et al.l l2009l) . scaling relations have 
to be established as a function of stellar pa rameters, 
such as luminositv ((Ch ristcnsc n-Dalsgaard fc Frauds^ 
1983; Kieldsen fc Beddinel ll995t iHoudek et al.l 119991 
Samadi et al. 2 0071 ). The influence of various parameters 
is tested, such as the evolutionary status and the stellar 
mass. 

Data, methods and definitions are presented in Sect. [5] 
and in the Appendix. After the comparison of the different 
methods used for measuring the global energetic parame- 
ters, scaling relations are derived in Sect. [3] The fine struc- 
ture of these relations is examined in Sect. 21 In Sect. [5l we 
propose a method for directly deriving the amplitudes of 
radial and non-radial modes, using the complete identifica- 
tion of the red giant oscillation spectra. This method is also 
used in Sect. |6] to compare observed bolometric amplitudes 
to available relationships, and in Sect.[7]to determine mode 
visibilities. 



2. Data and method 

2. 1 . Time series 
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Fig. 1. Local value of the exponent of the background 
as derived from Eq. |31 as a function of fmax- The dashed 
line indicates the median value estimated in the range [5, 
120 /iHz]. The increase of as at large frequency is an arte- 
fact due to the Nyquist frequency. The color code indicates 
the evolutionary status: clump stars in red, red giant branch 
stars in blue, and undetermined status in black. 

mainly due to large but rare interruptions. In a few cases, 
3-month long gaps occur for stars falling on the damaged 
CCD in the Kepler focal plane. 



2.2. Background and power excess 

To define the different components making up the signal, 
we used the classical phenomenological description of an 
oscillation power excess described by a mean Gaussian pro - 
file superimposed on the background (jMichel et al.ll2008h . 
This means that the power density spectrum is composed 
of three components: the background B, which is assumed 
to be dominated by stellar noise at low frequency, the os- 
cillation signal P centered at I'max, and white photon shot 
noise. The frequency z/max corresponds to the location of the 
maximum oscillation signal, according to the description of 
the Gaussian envelope: 



P{v) — -ffi/,„j,^ exp 



2a2 



(1) 



The st ars analyzed in this work have already been pre- 
sented ()Hekker et al.|[2011al and references therein). Their 

Kepler magnitudes are in the range 8 - 12.3, with most of ^yy JrLarvey-l iJte componen t! 
them fainter than magnitude 11. They show oscillations in Lorentzian (jHarvev^ll985^ : 
the frequency range [1, 280 /iHz], for which Kepler long- 
cadence data are well suited, having an average sample 
time of 1765.5 s, and so a Nyquist frequency of 283.5 /iHz. 
This work includes all available data, providing time se- 
ries of about 590 days for more than a thousand red gi- 
ants, giving a frequency resolution of about 19.7 nHz. The 
time series were prepared from the raw data (iJenkins et al 
I2OIOI ) according to the method presented bv iGarcia et al 



where 5varw = 2-\/21n2(T is the full- width at half-maximum 
(FWHM) of the mode envelope, and H^^_^^ the height at 
t'max- We note that, for the bright red giants considered 
here, the photon noise component is negligible compared 
to the other components. 

In the global approach, the background B is described 
by Harvey-l ike componen ts. Each component is a modified 



+ (27rt/ri)"' ' 



(2) 



()2011l ). which provides correction for various instrumental 
artefacts. The mean duty cycle is 91.5%, with data loss 



where Ti is the characteristic time scales of the i-th com- 
ponent of maximum height bi. The total number of com- 
ponents varies from 1 to 3 in the different methods. Such a 
model provides an accurate description of the background, 
although it is phenomenological. The introduction of a 
Lorentzian is linked to the assumption that the background 
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signal is composed of independent components with an ex- 
ponential decay in the time domain. In practice, the expo- 
ne nts aj, when fi x ed, ar e chosen to be close to 2 or 4. 

iMathur et alJ (j201lh explicitly addressed this descrip- 
tion of the red giant power spectra. Their work shows that 
the frequency (27rr)~^ of the granulation component close 
to fmax is, in fact, closely correlated with Vma-x', it varies as 
'^niax- III this work, we are mainly concerned with the value 
-^I'max of the background at J^max- In that respect, a local 
description of B{u) in the vicinity of t'max provides a model 
of sufficient precision. We have found that a polynomial 
approximation of the form 

BH^B^_(-^y" (3) 

is more precise than a linear fit. The expression holds 
in the frequency range [i/max - (^z^cnv, I'max + <5i^env]- The 
local parameters were estimated in the frequency ranges 
surrounding the region where oscillation power excess is 
observed. The exponent as is found to be — 2.1 ± 0.3 
(Fig, m . According to the relation between l/r and i^max 
(jMathur et al.|[201ll) . this supports an exponent close to 2 
in the Harvey profile (Eq. [2]), with a typical frequency l/r 
smaller than 2 fmax- If a modified expression including a 
term in i/"* is preferred to account for a rapid decrease of 
the background at higher frequency than z^max, the mea- 
surement of as close to —2 implies stringent conditions 
between Vmax and l/r. We conclude that, at lower frequen- 
cies, an exponent of 2 works reasonably well, whereas at 
higher frequencies a higher power of about 4 may be re- 
quired to reproduce a steeper gradient, presumably reflect- 
ing smaller length scales in the turbulent cascade describing 
the convection. 

The global description of the background (Eq. ^ is 
intended to be more accurate than the local description 
(Eq. in]). However, until we have a physical description of 
the background and its different scales, the global descrip- 
tion remains phenomenological. We must keep in mind that 
crosstalk between the different components of the signal 
might add or remove energy from the oscillation signal. 
In that respect, testing the local description of the back- 
ground, which is certainly not a correct solution, offers the 
possibility to fit the background with a minimum number 
of components. 

In order to obtain the global parameters of the Gaussian 
envelope, most of the methods used for this work first con- 
sider a smoothed spectrum. The best practice for achieving 
this step is described in Appendix and |B] 

3. Scaling relations for global energy parameters 

3.1. Comparison of the methods 

We focus on the three global parameters previously intro- 
duced for describing the power excess (Eqs. [TJ3]): 

— The FWHM diycnv of the Gaussian envelope, which mea- 
sures the extent of efficient excitation of the oscillations. 
As we already know from CoRoT observations, the 
FWHM increases significantly with i^max (Fig.H^). This 
envelope width is important to estimate the number 
nenv = Svonv/^i^ of radial orders with detectable am- 
plitudes. This number can be compared to the order of 
maximum signal, defined as ri„iax = i^max/A;^ — ei^Av). 



Table 1. Comparison of global seismic parameters 



method 


Af 






.x/10' 






A2Z 


0.269;0.755 


0.49;0.89 


0.90 


-2.10 


11.1 


-2.56 


CAN 


0.292;0.740 


0.71;0.82 


1.95 


-2.35 


5.75 


-2.31 


COR 


0.286;0.744 


0.74;0.84 


1.69 


-2.28 


3.76 


-2.23 


OCT 


0.271;0.756 


0.57;0.90 


1.74 


-2.38 


8.16 


-2.37 


SYD 


0.283;0.747 


0.63;0.85 


1.95 


-2.32 


5.36 


-2.39 



- Each column gives a doublet a; P, respectively coefficient and 
exponent of the power laws a v^^.^ of each parameter. 



- Frequencies are all expressed in ^Hz, and heights in 
ppm^/iHz~^. 

- The different methods provided results for a number of targets 
varying between 720 and 1220. 



— The oscillation height H is measured at fmax (Fig.[2jD). 

— The background signal underlying the oscillation signal 
is taken as the height Bi,^^^ at t-max (Fig. Ht). From 
B^ , we can derive the height-to-background ratio 

hZ1/Bu^.. (Fig.m). 

The preliminary normalization process (Appendix \K\ 
ensures all methods produced results in global agreement. 
Small differences nevertheless remain, mostly being re- 
lated to the way the background is modelled, particu- 
larly the number of components in the background pro- 
file. The differences have a stronger impact on the mul- 
tiplicative coefficients than on the exponents of the scal- 
ing relations. Despite the correction of the time-domain 
averaging effect towards high frequencies due to the half- 
hour integration time, expressed by the suppression fac- 
tor sine (7rfniax/2^'Nyquist), wc uotc that all methods are 
affected by the Nyquist frequency. When v^aax is close to 
I'Nyquist, it becomes difficult to estimate the background, 
and unsurprisingly all methods tend to underestimate it 
since fewer modes are observable. As a result, the following 
analysis is limited to stars with fmax < 200 /^Hz. 

The parameter showing the best agreement between the 
different methods is the envelope width (5z/env, while the pa- 
rameter with the largest spread is the height-to-background 
ratio. This refiects systematic variations in the way the 
different components of the Fourier spectrum were fitted. 
We note systematic differences in H^^^^ and B^^^^, often 
with anticorrelated variations. We also note a change at 
the clump-frequency, which is certainly real, and is studied 
in detail in Sect. HI Finally, we consider that the local de- 
scription of the background provides results intermediate 
between the different global descriptions based on Harvey 
profiles. We therefore chose to use it for analyzing the scal- 
ing relations in the following sections. 

3.2. Scaling relations 

Scaling relations for Suc-av, H^^^^^ and B^_^_^^ as a function 
of fmax show clcar global trends (Table [TJ Fig. Mean 
values of the coefficients and exponents of the best fits de- 
rived with all methods are given in Tabled The dispersion 
with respect to the power laws, even if larger than for the 
i^max — Av relation, are fairly small. Many reasons, with 
different origins, can explain this. Only physical reasons re- 
lated to the seismic properties are discussed below. Other 
factors, such as the infiuence of unidentified close compan- 
ions or background objects on the red giant light curves. 
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Fig. 2. Global seismic parameters as a function of t'max derived from the five methods used here, a) FWHM Svcnv of 
the Gaussian envelope; b) Mode height H,^^^^, as defined by Eq. [2 corrected for apodization; c) Background at i^max- 
The increase above 200 /iHz is an artefact due to the proximity of the Nyquist frequency, d) Height-to-background ratio, 
affected by the artefact on B,y^^^. The color code indicates the method used for the data analysis, as indicated in panel 
a). The deviations from the trends above 200yuHz are due to the proximity to the Nyquist frequency. 



are not relevant in this study since the red giants currently 
accessible were chosen to be bright and in uncrowded fields. 

— The scaling relation for Si/cnv observed for CoRoT tar- 
gets is confirmed. The 'saturation' above 200 /iHz is cer- 
tainly related to the difficulty in estimating the back- 
ground when i^max is close to i^Nyquist- As stated earlier, 
stars with t'max greater than 200 /iHz were not consid- 
ered in the following scaling relations. But we also note 
a saturation in the frequency range [100, 200 /iHz] that 
cannot be explained by this artefact. This effect is stud- 
ied in detail in Sect. SI 

— The number of observable degrees varies as (5fonv/Ai^. 
It dec reases slowly at lo w frequency since it scales as 
ly^-ll (|Mosser et al.l[20To h. 

— The stellar background signal varies as t'max''- The ex- 
ponent differs from the va lue —2 observed for subgiants 
and main-sequence stars (i Chaplin et aII[20TTh . 

— The exponents of the oscillation and background heights 
are approximately equal. As a consequence, the HBR 
appears to be constant with frequency (Fig. This 
was not the case in CoRoT data, where the background 
for many stars was dominated by the stellar photon 
noise, without a possibility to disentangle precisely the 
photon noise component from the stellar background. 



Here, photon noise is negligible compared to the stellar 
background for the majority of the targets. The HBR 
at fmax measures the energy density of the oscillations 
relative to the stellar background. The constant ratio, 
equal to about 3.8, is discussed in the next section. We 
note that the HBR values in Fig. [2ji have a high disper- 
sion in the frequency range [30, 60 /iHz]. The very large 
number of objects observed in this range excludes spu- 
rious variation due to a poor sampling: these stars def- 
initely present characteristics that cannot be described 
by a single power law in Vrnax- 

For the remainder of the discussion, all results are based, 
for consistency, on a single method (COR, Table ETTj) . 

4. Fine detail 

Since we lack theoretical explanations for the observed scal- 
ing relations, we examine in this section how these relations 
are inffuenced by stellar parameters. The role of the stellar 
mass is potentially im portant since the m ass distribution is 
degenerate with ^-max (jMosser et al.ll2010l ). So, we have de- 
rived an estimate of the asteroseismic mass from the scaling 
relations for Av and v^nax to address this mass influence. 

In a next step, we have examined whether the relations 
discussed in Sect. Ovary with the evolutionary status. The 
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Table 2. Scaling relations 



parameter 




coefficient a 


exponent /? 






0.276 ± 0.002 


0.751 ± 0.002 


^max — ^max 




3.52 ± 0.004 


0.253 ± 0.002 


SUcnv 




0.66 ±0.01 


0.88 ±0.01 






2.33 ±0.01 


0.13 ±0.01 






(2.03 ± 0.05) lO'^ 


-2.38 ±0.01 






(6.37 ±0.02) 10** 


-2.41 ± 0.01 



- Mean values from the combined results of all methods. 

- Each parameter is estimated as a power law of Vma.x, namely 
varying as az^maxi with a the coefficient and /3 the exponent. 

- Frequencies are all expressed in fi¥Lz, and heights in ppm^/^Hz. 

- Scaling relations concerning Sucnv are only valid for stars with 
I'max less than 100 /iHz. 
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Fig. 3. Comparison of two stars with similar large sepa- 
rations, one belonging to the RGB (KIC 4750456, Ai^ = 
5.89 /iHz, blue curve) and the other to the secondary clump 
(KIC 3758458, = 5.87 ^iHz, red curve). The x-axis is the 
reduced frequency n'; vertical lines indicate radial modes. 
The dashed curves represent 20 times the smoothed spec- 
tra; the regions corresponding to the FWHM are overplot- 
ted with solid lines. 



direct comparison of two stars with similar large separa- 
tions but different evolutionary status, one belonging to 
the red giant branch and the other to the red clump, clearly 
shows that the clump stars present much lower amplitudes, 
but larger values of t'max and Svcnv, hence a higher acoustic 
cutoff frequency (Fig. [3]). The relation between ^-max and 
the cutoff frequency has been shown observationally, and is 
now assessed theoretically (|Belkacem et al..,2011, ). 

4.1. Mass-radius relation 

In order to investigate how global seismic parameters vary 
with both mass and evolutionary status, we must illustrate 
how the mass, radius and evolutionary status are linked. 
Fig.|3]shows how the stellar masses and radii are distributed 
with J^niax and how they are correlated with the evolution- 
ary status. For stars on the red giant branch (RGB), the 
mass distribution is nearly uniform across the whole fre- 
quency range. On the other hand, the distribution of i^max 
and mass are clearly correlated for clump stars, as already 
shown with CoRoT data. The same features are clearly 
found in the mass-radius relation, since stellar radius and 
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Fig. 4. Asteroseismic mass as a function of i^max (panel a) 
and of the asteroseismic radius (panel b). The color code 
indicates the evolutionary status; clump stars in red, giant 
branch stars in blue, unknown status in dark grey. The pop- 
ulation of giants with low £ = I amplitude (identified later 
in Sect. 17.1.11) is indicated with black squares. The rectan- 
gles in the upper right corners indicate the mean value of 
the l-cr error bars. 



J^max are strongly correlated. A discussion in terms of stel- 
lar evolution, mass loss, and reset of the red giant structure 
at the heliu m flash for low- mass stars, has already been 
presented bv iMosser et al.l ([201l3). Here, the clear differ- 
ence between RGB and clump stars helps to identify the 
structure present in Fig. ^jp: the apparent slope change in 
the H^^^^ and B^^^^ relations is linked to the highly non- 
uniform mass distribution for fmax m the frequency range 
[30, 60 fiHz] due to clump stars. 

4.2. Power excess, mass, and evolutionary status 

We have revisited the scaling relations dealing with the 
power excess according to the evolutionary status of the 
red giants (Figs. OandEl). This reveals a clear difference 
between the two populations. 

— For stars with i^max m the range [25-50 /iHz]: red clump 
stars and RGB stars have similar masses, but the dis- 
tribution of the mass on the RGB is not correlated 
with fmax, whereas there is a clear correlation for red- 
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Fig. 5. Maximum height H^,^^^^ at i^max (panel a) and Fig. 6. FWHM in unit as a function of i^max (panel a) 

product H,y^_^_^^Si'cnv (panel b), as a function of t'max- and of the seismic mass (panel b). Evolutionary status and 

Evolutionary status and mean error bars are indicated as mean error bars are indicated as in Fig. |4l 
in Fig. m 



clump stars (Fig.|3]). At the same time, the global energy 
parameters are similar for both populations, but with 
somewhat different slopes (Fig. [SJ. 

— At higher i^max we find the secondary clump, consist - 
ing of He-burning stars of higher masses (lGirardil ll99i . 
With increasing mass, the maximum mode height of sec- 
ondary clump stars decreases significantly. In parallel, 
the FWHM of the Gaussian envelope increases. In stars 
burning helium in their core, the energy partition differs 
significantly: a wider range of modes is excited but with 
much lower oscillation heights. This can be understood 
when recalli ng the link between t ^ pax and the cutoff fre- 
quen cy Vc (|Kieldsen fc Beddind [TqqS iBelkacem et~all 
l201lh . At fixed Ai/ (as for the two red giants shown in 
Fig. 131) , a low-mass star has a much lower t'max than a 
member of the secondary clump. The explicit mass de- 
pendence of the FWHM Siyenv is shown in Fig. \E[ The 
product Hi.^^^Si'cnv (Fig. [5]d) tells us that oscillations 
in clump stars have less energy than in the RGB. 

— The saturation of dvcnv observed in the frequency range 
[100 - 200 /iHz] only concerns RGB stars. It occurs at 
too low a frequency to be related to the Nyquist fre- 
quency, and clearly departs from the scaling relation 
reported in Table [T] We have to conclude that the mech- 
anism of excitation is only efficient in a limited fre- 



quency range. This may be due to a change in the rela- 
tion between t^max and the atmospheric cutoff frequency 
(jBelkacem et al.ll201lD . Finally, we can derive a relation 
between Svonv and the stellar mass, expressed by the 
number of observable radial orders: 

^ Si^cnv _ / 2.72 + 0.87 M/Mq (clump) . . 
Av ~ \ 2-37 + 0.77 M/M© (RGB) 



4.3. Granulation background 

We have tested the new relation giving B^^^^ as 
a function of Ab =_(L/Lq)'^{M/Mq)-^{T^s/Tq)~^-^ 
(|Kieldsen fc Beddind[20Tlh . We made use of the asteroseis- 
mic estimate of the luminosity derived from the asteroseis - 
mic radius, the effective temperature (jBrown et al.l 120111) 
and the Boltzmann law. Such luminosity values certainly 
suffer from scaling uncertainties. This means that the ab- 
solute values are less reliable than relative variations. In 
other words, the exponent of the relation predicting B^^^^ 
is less biased than the multiplicative coefficient. 

The function B^_^__^^{v-at_s.yi) shows two different branches 
with different slopes (Fig. [7^), whereas ^^^^^^^(Ab) seems to 
be able to reconcile both evolutionary states (Fig. [7)d) . The 
phenomenological expectation seems verified, since the ex- 
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Fig. 7. Background -B;y„,,,^ as a function of fmax (panel a), 
and as a function of Bq Ab (panel b), with the solar value 
Bq — 0.19ppm^ /iHz~^. Evolutionary status and error bars 
are indicated as in Fig.^l The dashed line indicates the 1:1 
relation. 

oneii t is not far from unity, in agreement with lHuber et ahl 
201 Ih : 




(all) 

(clump) 

(RGB) 



(5) 



We note that the coefhcient of the fit of all red giants is 
very close to the value of the sol ar background at t^ max of 
about 0.19 ppm /iHz"i found bv lHuber et all (l20Tll) . This 
value has been used in Fig. [7)d. We see that the relation 
predicts too low values of the background at t'max- 

Having previously remarked that the HBR has a uni- 
form value across a wide frequency range, we have then 
to examine whether H^^,^^ also scales with Ab. This is the 
case, with exponents close to 1 for both RGB and clump 
stars, but with a different ratio. We measure: 



B, 



3.67 ±0.07 
4.02 ±0.12 



(clump) 
(RGB) 



(6) 



iKieldsen fc Beddind ()2011h postulated that the mean 
height of modes, observed in velocity, is proportional to the 
velocity power density of the granulation at i^max- Here, we 
show that this proportionality is verified in the photometric 



signal for re d giants. A similar s tudy for all stellar classes 
was made bv lHuber et al.l (|201ll) . who reached similar con- 
clusions. 



5. Amplitudes 



It has been 
tudes d epend 



suggested 

on the 

tio (IChristensen-Dalsgaard h Frandsen 

Kieldsen fc Beddind 119951: iHoudek et all 



that oscillation ampli- 
luminosity-to -mass 



1983 



1999; 



Samadi et al.l l2007l ). They also depend on the mode 



degree. As we found a large diversity in the red giant 
oscillation spectra, we propose in this section a method for 
measuring separately the radial and non-radial amplitudes 
in red giants facilited by the automated mode identification 
provided by the red giant universal oscillation pattern 
(Appendix jCj. The measurements will then allow us to 
derive bolometric amplitudes and also mode visibilities, in 
Sect. |S| and |7| respectively. 

5.1. Radial amplitudes 

Previously, the amplitudes of radial modes have been de- 
rived from the mean amplitude in the oscillation enve- 
l ope divided by the t otal visibility factors of the modes 
()Kieldsen et alj 120081: iMichel et aLll2009D . Here, we have 
computed spectra as a function of the reduced frequency n' 
(Eq. 



IC.2|) . Squared amplitudes are then given by integrat- 
ing the height across the mode width, as estimated from 
the fit of the ridges, with the integral being corrected for 
the background contribution. That is, we have computed 
the squared amplitudes for each radial order from: 



Al(n) 



5v 



\p{n) — i?] An 



(7) 



with p the power density spectrum, 5v the frequency res- 
olution, and B the local background. The boundaries 620 
and eo3 bracket the frequency range where radial modes are 
expected: the radial mode (n, 0) lies between the {n — 1, 2) 
and (n — 1, 3) modes. They have been defined according to 
the parametrization proposed bv lMosser et al.l pOllb ). We 
have checked that the measures of the radial amplitudes are 
stable with respect to slight modifications of these bound- 
aries. The major contribution to the uncertainties comes 
from uncertainty on the background correction, resulting 
in a mean precision of about 10 %. 

5.2. Non-radial amplitudes 

Non-radial amplitudes have been calculated in the same 
way, with the appropriate boundaries eui between the ad- 
jacent degrees i and i': 



A]{n) 



n-\-e^^/i 



[p{n') - B] dn' 



(8) 



We have used the mean values ei2 — —0.22, 620 — —0.065, 
eo3 = 0.17, and 631 = 0.27 derived froni the u niversal 
red giant oscillation pattern ([Mosser et al.ll20lTbl) . These 
boundaries are globally shifted by 0.008 (n — nmax) to ac- 
count for the mean curvature of the spectra, and modulated 
in large separation according to the ex act location of the 
ridges (Table 1 of iMosser et all 1201 Ibf ). In case of mixed 
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Fig. 8. Identification of tire ridges in red giants observed by Kepler. Targets were cfiosen with oscillation spectra exhibiting 
very different visibility coefficients (see Sect. [7]). Left side: high visibility of ^=0 to 2 from top to bottom. Right side: low 
visibility of i—0 to 2 from top to bottom. The spectra are corrected for the background contribution (negative values 
being omitted), and plotted as a function of the reduced frequency. The color code, derived from Eqs. [7] and[51 presents 
radial modes in red, = 1 in blue, £ = 2 in green, and = 3 in light blue. 



modes, the amplitude of a certain degree we measure cor- 
responds to the sum of all individual components of a given 
pressure radial order n. A few examples are given in Fig. [S] 
From the definition given by Eq. [51 it is clear that, 
in case of rotational multiplets, the method integrates the 
squared amplitudes of all components. In other words, the 
method is not affected by the unknown stellar inclination. 

5.3. Mean values of individual amplitudes 

The mean value of the individual amplitudes has been ob- 
tained under the assumption that the energy partition is 
Gaussian, following Eq. [TJ 



J2 ^'(^) / J2 



2ct2 



(9) 



Formally, one would expect the weighting to be performed 
before the summation, but such a method is too sensitive 
to the varying amplitudes created by the stochastic excita- 
tion. Therefore, with Eq. [SJ we intend to estimate a mean 
value not affected by this source of noise. Considering five 
radial orders is in agreement with the observed width of the 
envelope and allows us to pick the highest peaks. 



5.4. Individual and globally averaged radial amplitudes 

The radial amplitude has also been directly determined, by 



(10) 



We have compared (Aq)^/^ to the radial amplitude Aq, 
and shown that both amplitudes are closely linked, with 
Aq ~ 1.07 (Aq)^/2. However, to avoid scattering due to 
the finite mode lifetimes and for coherence of the visibility 
calculation, we have decided to consider the radial ampli- 
tude (Aq)^/2. This value can be compared to the non-radial 
amplitudes in order to compute the mode visibilities. It 
can also be compared to the global a verage amplitudes de- 
rived from the tot al mode visibility ([Kieldsen et al.l [20081 : 
iMichel et alll2009h 



vi 



(11) 



The examination of the total mode visibility Vj-^j in Sect. [7] 
allows to compare {A) and [A^Y^"^ in more detail. 

6. Bolometric amplitudes of radial modes 

We now use the bolometric co r rectio n performed for Kepler 
observations by iBallot et al.l (|2011f ) to translate the ob- 
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Fig. 9. Bolometric amplitude Aho\, in Ppm, as a function of 
(panel a), Aq [LMq/ L^Mf-^T^s /Tq)-''-^ (panel b) 
and Ag (i/LQ)(M/M0)-i-^(T/r0)O-5(Teff/To)-2-75 (panel 

c), with the solar value Aq ~ 2.53 ppm. Evolutionary sta- 
tus and error bars are indicated as in Fig. The dashed 
lines show the 1:1 correspondence. 



served radial amplitudes to bolometric amplitudes. We then 
examine how the bolometric radial amplitudes scale with 
various parameters and test the available theoretical am- 
plitude scaling relations. 



6.1. Bolometric correction 

The bolometric correction provided by iBallot et al.l (|2011f) 
allows us to derive bolometric amplitudes from the radial 
amplitudes: 

Aoi = {AlY'^ Cboi with Cboi = (T;ff/rK)"-«° (12) 

with Tk — 5 934K. This correction accounts for the wave- 
length dependence of the photometric variation integrated 
over the Kepler bandpass. The variation of Aho\ with i/max 
shows clearly the differences between RGB and clump stars 
(Fig. [Sh-) ■ We note that the difference is much more pro- 
nounced for the secondary clump, which contains objects 
with higher mass than the primary clump. 

6.2. Comparison witli (L/Af)"'^ 

We have compared these amplitudes to different relations 
depending on the stellar luminosity and mass. As in Sect. 01 
both L and M were derived from the usual asteroseismic 
scalings relations. As discussed above, such scalings avoid 
the use of grid modelling, which only works at fixed physics, 
but certainly lower the quantitative validity of the outputs 
since the scaling relations are not yet accurately calibrated. 

We have t e sted t he variation in (L/M)" '' proposed by 
iSa madi et al.l (|2007[ ) for Doppler measurements, and ver- 
ified that the observed exponent is effectively near the 
val ue 0.7 derived from 3-D simulations. As alre ady shown 
bv i Mosser eraLl ()2010[ ) and iHuber et al.l (|2010t ). we find a 
shghtly higher exponent, of about 0.8 (Fig.jHb and Table|3]). 
We have introduced the correction from Doppler to photo- 
metric measures, expressed by a factor \/Tcs, under the 
assumption that the oscil lations propagate adiabatically 
(jKieldsen fc Beddind 119951 ). We also adopted in the scal- 
ing relation the bolometric amplitude Aq ~ 2.53 ppm of 
sol ar radial modes observed in broad-band photometry 
(Mi chel et al.l [20091) . We note that the relation underesti- 
mates the amplitude by a factor of about 2, and that the 
exponents are different for the different evolutionary states 

(Fig.ib). 

Th e theoretical relation proposed by ISamadi et al.l 
(pool was derived for main-sequence stars, and applying 
it to red giants is a large extrapolation. Indeed, investi- 
gating mode driving in red giants would require a non- 
adia batic treatment. F rom the oscillation energy equation 
(e.s. lUnno et al.lll989L Chap. IV, Eq. 21.14), dimensional 
arguments show that non-adiabatic effects scale as the ra- 
tio L/M, which is roughly 50 times greater for red giants 
than for main-sequence stars. We might therefore expect 
the physics underlying mode driving to be different for the 
two cases. In principle, non-adiabatic effects must also be 
considered when converting between velocity amplitudes 
(as provided by theoretical computations) and photometric 
amplitudes (as measured with Kepler). However, in the ab- 
sence of a reliable non-adiabati c theory, we have to adopt 
the adiabatic conversion (e.g. iKieldsen fc Beddind 119951: 
ISamadi et al.ll2010l) . These points deserve thorough theo- 
retical investigation, but this is beyond the scope of this 
paper. 

6.3. Comparison with LjM^-^ 

We also t ested the revised amplitude scaling relation sug- 
gested bv IKieldsen fc Beddind ()2011|) . It includes not only 
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Table 3. Bolometric 


amplitudes 




2.0 


- 


scaling with 


coefficient 


exponent 






i^max all 


860 ± 30 


-0.71 ±0.01 


1.5 


- 


clump 


4700 ± 500 


-1.18 ±0.03 






RGB 


650 ± 60 


-0.63 ±0.02 






L/M all 


3.5 ±0.2 


0.79 ±0.01 






clump 


0.92 ±0.11 


1.12 ± 0.04 


^" 1.0 




RGB 


4.8 ± 0.4 


0.71 ±0.02 






Aa all 


8.2 ± 0.2 


0.70 ± 0.01 


< 




clump 


5.3 ±0.3 


0.83 ± 0.02 






RGB 


9.8 ±0.4 


0.65 ±0.02 


0.5 




the mass and luminosity, but also the effective temperature 






and the mode lifetimes r: 




0.0 





•+,•*• .** V 




Aho\ oc Aa with 



(L/io) (r/ro) 



0.5 



iKieldsen fc Beddind (1201X11 . As a consequence, this theoret- 
ical relation largely overestimates the observed bolometric 
radial ampli t udes, as also shown bv lHuber et al.l ()201lD and 
IStello et all (1201 Ih . This means that the proposed relation 
does not provide a complete physical explanation. We also 
note that the prediction of ^boi proportional to L is ob- 
servationally unlikely. For red giants, L is the most rapidly 
varying parameter in Eq. [T31 so that the observed slope 
in Fig. almost directly translates into the luminosity 
exponent. The agreement between observed and predicted 
values would require r to be approxi mately proportional 
to Av, which is cle arly not observed (jHuber et al.l l2010t 
iBaudin et a^l2011a^ . 

Despite the quantitative disagreement with theoretical 
predictions, we see from Figs. [Sb and c that an empiri- 
cal scaling relation can describe the variation of Aho\ as a 
function of the stellar luminosity, mass and effective tem- 
perature. Such a fit based on red giants in open clusters, 
hence with a n acc urate mass estimation, was proposed by 
IStello etUri (|2011h . which agrees with a fit based on red 
giants, subg i ants a nd main-sequence stars performed by 
iHuber et all (|2011h . 



-1.0 



-0.5 



0.0 



0.5 



1.0 



(M/M0)i-5 (Toff/T0)2-75- (13) 

iBaudin et al.l ()2011aD explored the observed mode lifetime 
for CoRoT data. They suggested that for red giants, r is 
about 30 days and varies as T~^'^^^'^. We use these esti- 
mates here. When expressed as a function of Aa, the two 
branches corresponding to RGB and clump stars show a 
closer correlation (Fig. Et), so that the distribution of ^boi 
as a function of Aa is less dispersed than the distribution as 
a function of (L/M)°-^. However, the best fit gives a vari- 
ation Aho\ cx a" ''"^^" ' "^ ins tead of ^boi oc Aa proposed by 



Fig. 10. Variation in metallicity of the bolometric ampli- 
tude scaled to a best fit in L/AI^-^. Same color code as in 
Fig. m The red and blue lines corresponds to the linear fits 
i n Z. T he uncertainties on Z are obtained from lBrown et all 
(120 llh . 
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Fig. 11. Visibility Vi as a function of I'max, with the same 
color code as in Fig. |4l Large black symbols indicate the 
population of stars with very low values. 



7. Spatial response functions 

We have examined the visibility of each degree, obtained 
from 



6.4. Influence of metallicity 

It has been suggested t hat oscillation ampl i tudes depend 
on t he stellar metallicity (|Houdek et al.|[T999t ISamadi et al.l 
I201Q ). In order to test this dependence, we have plotted 
the bolometric amplitudes as a function of the metallicity 
(Fig. nop . This comparison requires normalized ampl itudes, 
which were obtained using a fit in L/M^-^ (Hube r et al.l 
1201 Ih . The met allicity value s were taken from the Kepler 
Input Catalog (jBrown et al.ll201l[ ) and have uncertainties 
of about 0.4 dex. We found that an increase of 1 dex in 
Z gives a moderate amplitude increase of about 10±5%. 
Further work, based on improved values of the metallicity, 
will be necessary to make the link more precise. 



{AD I {AD, 



(14) 



with (Aq) given by Eqs. [7] and [SI and obtained in 

a similar way. For pure p modes in a 4800 K red giant, 
one expects spatial response functions in power to vary 
as ~ 1.5 4, V:^ - 0.58 and V^} ~ 0.043. These values, 
derived from lBallot et al.l (|201ll ). assume energy equiparti- 
tion between the different degrees and take into account the 
influence of the limb darkening. They are expected to de- 
crease when the effective temperature increases (Table |4]). 
Fig. [5] shows the background-corrected spectra of different 
red giants exhibiting different mode visibilities. 
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Fig. 12. Visibility of the £ = 1,2 and 3 modes (from top to 
bottom), as a function of the effective temperature. Same 
color code as in Fig. [TTl Dashed lines are the fits derived 
from I Ballot et al.l ()2011.) . Blue and red solid lines are the 
fits for RGB and clump stars, respectively. 



7.1. Dipole modes 

The visibilities Vi of the dipole modes are plotted in Fig.fTD 
as a function of t'max and in Fig. 112b . as a function of the 
effective temperature. We first note a large dispersion of 
the values and a discrepancy between observed and pre- 
dicted values. In general, the observed Vi appear to be 
smaller than expected. Despite including the mixed modes 



Fig. 13. Visibility of the i = 1 and 2 modes, as a function 
of the seismic mass. Same color code as in Fig. [TTJ 



in a very broad frequency range, in fact larger than Av/2, 
we observed (V^^) = 1.34 ± 0.02 instead of the expected 
value of 1.54 (Tabic |4]). Individual uncertainties on are 
about ±0.08. They are principally due to the background 
modelling and to the difficulty of automatically identifying 
all mixed modes. 



7.1.1. Very weak Vi values 

We note the presence of a family of stars with a very weak, 
almost absent 1=1 oscillation pattern. It mostly comprises 
less evolved stars, with z^max > 50 /iHz and < 0.8, which 
appear clearly as outliers. As a consequence of the low £ — 1 
amplitudes, their evolution status is not determined, except 
for four stars, two in the RGB and two in the clump. Very 
low Vi values are also measured at the clump, but without 
the clear difference compared to other stars. Therefore, we 
do not include them among the set of stars with very weak 
Vi values. We limit it to fmax > 50 /iHz and a low Vi 
(about 12.5 ^max^^j with :^max in /iHz). 

Different tests have shown that, apart from a low £ ~ 1 
visibility, these stars mostly share the common characteris- 
tics of other stars. Unsurprisingly, they have lower heights 
than expected from the empirical scaling relation since 
£ = 1 modes are almost absent, but their bolometric am- 
plitudes are close to the mean values. More precisely, if one 
believes that the location in the y^boi(i'max) relation (Fig.lH]) 
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makes it possible to identify the evolutionary status, clump 
stars with a low Vi show normal bolometric amplitudes, 
but RGB stars with a low have lower amplitudes than 
expected. This could be related to the fact that these stars 
have slightly higher masses compared to the mean values 
(Fig. ll3b .). However, some stars do present simultaneously a 
high asteroseismic mass and a normal Vi . Their £ — 2 visi- 
bilities are slightly lower than average, but not significantly 
so. Furthermore, these stars are uniformly distributed in 
temperature and metallicity. Therefore, we have to con- 
clude that only i = 1 modes are affected. One may imagine 
that their suppression results from a very efficient coupling 
between pressure and gravity waves. This efficient coupling 
yields a very high mode mass, hence a very low observed 
amplitude at the surface. Examining the causes of very low 
Vi values will need further work. It may also help to un- 
derstand the non-identification of £ = 1 modes in previ- 
ous observations o f high-mass giants ()Frandsen et al.ll2002t 
iBaudin et alllMlb) . 

7.1.2. A surprising energy equipartition 

Having identified the population with abnormally low Vi 
values allows us to calculate the mean £ = 1 visibil- 
ity of 'normal' red giants: it remains lower than the ex- 
pected value (1.38, versus 1.54), and uncertainties are un- 
able to explain the discrepancy. We note that the lower- 
than-expected values are independent of fmax, hence 
independent of the evolution. Coming back to the way the 
amplitudes were determined, we see that the sum of the 
squared amplitudes of all = 1 mixed modes correspond- 
ing to a fixed radial order is only slightly less than the total 
expected squared amplitude. We show later that the deficit 
can be explained by the energy being spread far away from 
the expected location of the £ — 1 p modes, which are 
not being detected by the automated method presented in 
Sect. El 

In terms of coupling, the energy injected in the pressure 
waves near the stellar surface seems to be shared among all 
mixed modes associated to a given radial order. The en- 
ergy equipartition seems valid when all squared amplitudes 
of the mixed modes associated to a given radial order are 
summed. Therefore, in terms of mode mass, one derives the 
equivalence between the radial and non-radial modes: 

1^1 

with rig being the gravity order of the g modes mixed with 
the p mode of radial order n. Each mixed-mode mass Mn^ 
is much higher than the radial mode mass, as discussed by 
iDupret et al.l (|2009t ). but the total as defined by Eg. [TSlis 
close to the radial mode mass. Future work will determine 
whether this is a coincidence or based on physical princi- 
ples. 

7.1.3. Variation of with Tcs, M and metallicity 

The variation of Vi with temperature matches the expec- 
tation, although with large uncertainties (Fig. [T^ . Table 
H]). The visibility Vi is supposed to increase towards low 
temperatures. If the extrapolation is valid, one may expect 
the observed £ — 1 modes to be more dominant for red 
giants towards the tip of the red giant branch. For these 



Table 4. Mode visibility 



visibility 


population 


value at 4 oUU K 


10 dVi /dicff 




expected 


1.04 


— U.UO 




clump 


1.34±0.02 


-0.07 ±0.09 




RGB 


1.35±0.04 


-0.10 ±0.12 




expected 


0.58 


-0.07 




clump 


0.59±0.01 


-0.02±0.06 




RGB 


0.64±0.02 


-0.03±0.06 


vi 


expected 


0.036 


-0.02 




clump 


0.056±0.01 


±0.03±0.01 




RGB 


0.071±0.01 


±0.04±0.01 




expected 


3.16 


-0.15 


clump 


2.98±0.02 


-0.07±0.10 




RGB 


3.08±0.04 


-0.11±0.15 



stars, the oscillation spectra cannot yet be analyzed since 
the frequency resolution is not high enough. 

We have also examined how the spatial responses 
vary with the stellar mass (Fig. [T3k). We remark that the 
high values of are concentrated at intermediate mass, 
between 1 and 1.7 times the solar mass. The visibility of 
£ = I modes is, on average, smaller for higher mass. The 
variation of Vi with metallicity shows no significant cor- 
relation, but this result may be due to the very uncertain 
metallicity determinations in the Kepler Input Catalog. 

Since correlations between Vi and any fundamental 
stellar parameters or seismic global parameters just dis- 
cussed are not very strong, we suspect that the scatter in 
Vi is related to the conditions that govern the coupling of 
pressure and gravity waves that contribute to the mixed 
modes. These conditions appear to be highly sensitive to 
the locations of the inner and outer cavities where g and p 
waves propagate, respectively. Future analysis will have to 
show if one can use the differences in visibility by correlat- 
ing them with properties of the eigenfrequency spectrum 
and drawing conclusions about the interior structure. 

7.2. Quadrupole modes 

The visibilities V2 of the quadrupole modes are plotted in 
Fig- nib- Compared to dipole modes, the dispersion is much 
lower. Uncertainties on V2 are smaller, about ±0.05, since 
the integration interval is reduced, and since uncertainties 
in the background modelling have less influence. 

The mean value is very close to the expected value. This 
is surprising, since £ — 2 modes are also mixed and so, as for 
dipole modes, one would expect lower V2 values. However, 
with a narrower g-mode spacing than for £ — 1 modes, in 
agreement with the asymptotic expectations, the density of 
£ = 2 mixed modes is much larger than for £ = 1. Perhaps 
this larger density, or a better coupling, compensates for 
the larger mode mass. This agrees with the better trapping 
of quadrupole modes, compared to dipole modes, expected 
by IDupret et al. (2009). 

We have checked that the variation of V2 with the ef- 
fective temperature also agrees with the theoretical expec- 
tation. Contrary to V^, the mass dependence of V2 seems 
to be nearly flat (Fig.flSb). 
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Fig. 14. Visibility V^^^. as a function of J^max- Same color 
code as in Fig.[TT] The fits of RGB, clump and low Vi stars, 
in blue and red solid lines or black dashed line, respectively, 
are computed for t'max ^ 50 /iHz. 
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Fig. 15. Ratio (A)/(Aq)^/^ as a function of J^max- Same 
color code as in Fig. [TT] The fits of RGB, clump and low 
Vi stars, in blue, red, and black, respectively, are valid at 
i^max > 50/zHz. The dotted line indicates the fit at fre- 
quency lower than 30/zHz. 



7.3. £ = 3 modes 

The visibilities oi £ — 3 modes have been calculated in a 
similar way (Fig. [Et). Absolute uncertainties on are 
smaller than for and V2, about ±0.03, but relative un- 
certainties are much higher due to their low intrinsic ampli- 
tudes. Individual checks confirm the reality of £ = 3 peaks 
with large amplitudes. However, we also identify the spuri- 
ous contribution of g-dominated £ = 1 mixed modes, which 
are located far away from the expected position of pure 
i = I pressure modes. They can contribute a significant 
part of the energy automatically assigned to the £ — 3 
modes. This may explain part of the deficit of the £ — 1 
energy reported above. 

We note that the vast majority oi £ — 3 modes are ob- 
served as single narrow peaks. This should indicate that, 
most of the time, only one £ = 3 mixed mode per Av- 
wide frequency interval is efficiently excited. Further obser- 
vations with a higher frequency resolution will help to make 
this point more clear. 

7.4. Total visibility and radial amplitudes 

We have estimated the total visibility Vj-Q^ = J2'^=o ^"^^ 
note a large spread of the values, with RGB and clump 
stars having different mean values (Fig. [H]) . We also note 
that the star s wit h a low also show a low total visi- 
bility. [Ballot et al.l (|201l[ ) have calculated red giant limb- 
darkening functions, and derived a mean value Vtot - 3.18 
for Teff = 4 800K and solar metallicity. We find a mean 
value of about 3. 06, in better a greement with the value 
3.04 calculated bv [Bedding et all l}2010). 

The mean value of V^^^. can be used to improve the coeffi- 
cient defining the globally averaged amplitude (^) given by 
Eq.[TlJ However, this globally averaged value is only indica- 
tive, due to the dispersion in V^^^ fFig. [Ti)) . This limits the 
use of {A) for deriving precise red giant radial amplitudes. 
Therefore, we discourage using the globally averaged value 
(A) for estimating the bolometric amplitude of a given red 
giant. Using (A) instead of (Aq)^/^ for establishing scal- 
ing relations is possible, but with an uncertainty on the 



exponent in L of about 0.04. This value results from the 
differences in the slopes observed in the scaling relation of 
(A)/(ylQ)^/^ to z^max (Fig.[in]), taking into account the fact 

that L OC fmax- 

8. Conclusion 

With Kepler photometric data recorded up to Q7, provid- 
ing 590-day long time series, we have determined scaling 
relations of the global parameters that describe the oscilla- 
tion power excess of red giants. About 1200 red giants were 
analyzed. 

8.1. Global parameters and scaling relations 

We have compared different methods, and shown that the 
background modelling infiuences the results. This infiuence 
does not prevent firm conclusions, but certainly needs fur- 
ther work to fully disentangle the fit of the granulation 
background from the oscillation power excess. There are 
indications that Harvey components with an exponent of 
2 are sufficient for modelling the background of red giant 
spectra in the frequency range around t^max- 

We have found a significant influence of the red gi- 
ant evolutionary status on the scaling relations describ- 
ing the as teroseismic and funda mental stellar parameters. 
Following [Mosser et al.l (|2011a[ ). we show that red-clump 
stars have a mass distribution clearly linked to the seismic 
parameter ^'max, hence to the stellar radius. Therefore, the 
scaling relations expressed as a single power law of i^max 
present a large dispersion around the clump. 

As a consequence of the different mass distributions, red 
giants on the RGB have larger oscillation amplitudes than 
the clump stars. Members of the secondary clump have dif- 
ferent properties: compared to RGB stars, they have much 
higher masses, lower oscillation amplitudes, but with more 
modes excited. 

We have compared the energy in the oscillations to its 
counterpart in the stellar background. The mean height in 
the oscillations is a fixed fraction of the energy density in 
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the background at fmax- This fraction depends only slightly 
on the evolutionary status. Since both granulation and os- 
cillation are due to convection, such a result clearly indi- 
cates that the excitation mechanism of the oscillation has 
the same efficiency for a very large variety of red giants. 

For both the background height and the bolometric am- 
plitude, the agreement between the RGB and clump stars 
is better when considering in the scaling relations a mass 
dependence not exactly opposite to the luminosity depen- 
dence. The bac kground height seems t o scal e as L^/M^T^g, 
as predicted bv lKieldsen &: Beddind (j201lt) . 



8.2. Bolometric amplitudes and mode visibility 

We have proposed a new method for deriving radial and 
non-radial amplitudes, and we have shown that the ampli- 
tudes obtained with a global integration are not appropriate 
for red giants, due to the presence of mixed modes that can 
significantly modify the total mode visibility. 

Radial amplitudes have been translated into bolometric 
amplitudes for comparison with theoretical expectations. 
We have shown that current theoretical scaling relations are 
unable to provide an acceptable fit to the observed photo- 
metric amplitudes. There are strong observational indica- 
tions in favor of an exponent in luminosity in the range [0.7, 
0.8]. A negative exponent of the stellar mass larger, in abso- 
lute value, than the luminosity exponent helps to reduce the 
discrepancy between R GB and clump stars ()Huber et al.l 
[20TltlStello eraII[20TTh . but it seems that ingredients other 
than L and M are necessary to properly account for the ob- 
servations. 

We have also derived the mode visibility of £=1, 2 and 3 
modes. The presence of mixed-modes reduces the observed 
values compared to expectations. As a result, we observe 
a new type of energy equipartition, where the squared am- 
plitude of the pure pressure modes (as radial modes) seems 
to be shared by all mixed modes. We note a high disper- 
sion that cannot be related to global parameters, and have 
identified a class of objects with very low £ = 1 amplitudes. 
It probably comprises red giants more massive than 1.3 Mq 
where the coupling between pressure and gravity waves is 
so efficient that all £ = 1 mixed modes have a very high 
mode mass. 



Appendix A: Normalization of the power density 
spectra 

Different methods have been used an d compared 
for t he analysis of the p o wer e xcess (iHuber et al.l 
20091: iMosser fc Appourchauxl [20091: iHekker et al.l l2010t 



Kallinger et al.l 120101: iMathur et al.M2010[ ). A comprehen- 
sive comparison of th e properties and chara cteristics of the 
methods was given bv IHekker et al.l ()2011aD . but only deah 
with the large separation and the frequency f,nax of 
maximum oscillation signal. A comparison of complemen- 
tary an alysis method s appli ed to the Kepler short-cadence 
data bv lVerner et al.l ()2011t ) presented results on the max- 
imum mode amplitude, but not for red giants. For the cur- 
rent work on red giants, parameters characterizing the os- 
cillation power excess have been compared. 

To ensure a correct comparison of results obtained with 
different methods, it was first necessary to normalize the 
outputs. The computation of power density spectra per- 




35 40 
(/^Hz) 

Fig. A.l. Infiuence of the width of the smoothing param- 
eter for a typical red-clump giant (KIC 1161618). The dif- 
ferent colors indicate the different values of the smoothing 
(in units of Av). 

Table A. 2. Influence of the smoothing on the seismic 
global parameters of KIC 1161618. 
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1.00 
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0.99 


1.00 
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1.28 


HBR(5)/HBR(1) 


1.01 


1.00 


0.81 


0.73 


0.58 


as(<S) / aB(l) 


0.99 


1.00 


1.10 


1.13 


1.16 



formed with Lomb-Scargle periodograms has taken into ac- 
count the correction for the duty cycle of each target, in 
order to obtain spectra with the frequency resolution corre- 
sponding to the total observation time and frequency up to 
the Nyquist frequency (i^Nyquist = 283.2 /iHz) related to the 
mean time sampling i5tLc of the Kepler long-cadence data. 
We chose to compute power density spectra (PDS) with 
the following normalization: a white noise signal recorded 
at the Kepler long-cadence sampling (5tLc with a noise level 
at (in ppm) gives a PDS with a spectral density a^, such 
that = 2a^5ti^c/N where N is the number of points in 
the time series. The characteristics of the methods used for 
this work are briefly described in Table lA.ll 



Appendix B: Smoothing 

Most of the methods use a smoothed spectrum to obtain 
the global parameters of the Gaussian envelope (Sect. [3]). 
Since the width of the G aussian envelope of red giant oscil- 
lation (Eq. ID) is narrow osser et al]|2010l) , this step must 
be performed carefully. An example for a typical red-clump 

and 



star is given in Fig. lA.ll The parameters H^^ 
have been calculated for the local description of the back- 
ground smoothed with different filter sizes and are summa- 
rized in Table IA.2I We have used a Gaussian window for 
performing the smoothing and have tested different values 
S of the FWHM. The mean value of the large separation 
acts as a natural characteristic frequency, so we expressed 
the FWHM in units of Ai/ and have tested values vary- 
ing from 0.75 Ai/ to 2.0 At^. When this value increases, the 
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Table A.l. Characteristics of the different methods 



pipeline A2Z° CAN" COR'^ OCT" SYD'= 

spectrum Lomb-Scargle periodogram — > power density spectrum (PDS) 

frequency axis corrected for the duty cycle 

frequency sampling Su — l/T, with T — N Sti^c the total observation duration 

i^yquist ~ 0.5/(5tLC, with Sti^c the Kepler long-cadence sampling 

PDS normalization =2 at 5ti,c/N 

i^max centroid Gaussian 

or peak of the 
smoothed spectrum 

smoothing 1 — Az/ Gaussian none 1 — Az^ Gaussian filter 

fihcr 

background global global local local global 

Harvey Harvey oc second order Harvey modified 

polynomial 

value of exponent free 4 free 2 2 and 4 



number of components 2 3 112 

References: a) iMathur et all ||2010D : V) iKallinger et al.l (|20ia ): c) iMosser fc Appourchauxl ||2009|): d ) iHekker et all l|2010l ): e) 
iHuber et al.l (|2009l 'l. For more information on the modelling of the background. see lMathur et alJ (|201ll ^. 



measured value of t'max decreases, with variations as large 
as 4%, depending on the method. In parallel, Hi,^^^ de- 
creases and B,y^^^ increases when the filter width increases, 
which results in variations of the height-to-background ra- 
tio (HBR) larger than a factor of two. As expected, Jfcnv 
also increases with S. This variation can be easily under- 
stood by recognizing that increasing the smoothing will 
spread the power of the oscillations relative to the back- 
ground. We consider an acceptable compromise to be a 
Gaussian filter with a FWHM equal to the mean large sepa- 
ration Av. Such a width is large enough for smoothing the 
influence of individual contributions of different degrees, 
and narrow enough to avoid the dilution of the global pa- 
rameters, as shown by the example given in Table [5^2] when 
the parameter S gets larger than 1, all terms show signifi- 
cant variations. 

This study shows the importance of the smoothing pa- 
rameter. A systematic variation induced by the smoothing 
as large as 4 % in ^'max is important because this parameter 
plays a key role in sca l ing relations for s eismic parame- 
ters (jMosser et al.l 120101: iHuber et aLll2010D . Since i/^ax is 
used f or deriving the seismic mass and radius (iMosser et al.l 
[20Tot IKallinger et aD l20Tot IHekker et al.1 |2011bO, a bi^ 
should be avoided in order to avoid subsequent biases in 
the stellar parameters. A systematic relative error cr^^^^ on 
f^max translates into a bias in the determination of the de- 
rived stellar parameters, of the order of (Ti,^^^ and Scr^^^^^^ 
for, respectively, the relative precision of the radius and 
mass. 

A variation of a factor of two of the HBR depending 
on the smoothing means that studying the partition of en- 
ergy between the oscillation and the background requires a 
careful description. According to the solar example, longer 
time-series are certainly required to reach the necessary pre- 
cision to draw firm conclusions. 



Appendix C: Mode identification 

The determination of the mode visibility (Sect. [J] and [7]) re- 
quires the complete identification of the p-mode spectrum. 
This is automatically given b y the r ed giant universal oscil- 
lation pattern (Mosser et al. 2 011bD . with the parametriza- 
tion of the dimensionless factor s of the asymptotic de- 



velopment ()Tassoull ll980'). The function e{Ai/) enables the 
identification of the radial modes 

iyn,e=o = [n + e{Ai^)] Ai^, (C.l) 

with Ah' the large separation averaged in the frequency 
range [i^max - '^I'onv, ^'max + Sugnv]- For simpHcity, we in- 
troduce the reduced frequency, dimensionless and corrected 
for the s term of the Tassoul equation: 

n' = v/Aiy-e{Av). (C.2) 

According to lHuber et all (|2010 D and' Mosser et al.l (|2011bD . 
e is mainly a function of the large separation. For radial 
modes, the reduced frequency is very close to the radial 
order, except for possible very small secondary-order terms 
in e{Ai') that do not depend directly on Ai/ and do not 
hamper the mode identification. 

The determination of the evolutionary status of the gi- 
ants (Sect. 14. 2[) . namely the measurement of the g-mode 
spacing of £ = 1 mixed modes, is then achieved wi th the 
automated method described bv lMosser et al.l (|2011al) . This 
could be done for about 674 targets out of 1043 with high 
signal-to-noise ratio time series. This method measures the 
g-mode spacing of the mixed modes from the Fourier spec- 
trum of the oscillation spectrum windowed with narrow fil- 
ters centered on the expected locations of each pure £ = 1 
pressure mode. 
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